arXiv:1503.01775v2 [astro-ph.CO] 7 Apr 2016 


The EFT of Large Scale Structures at All Redshifts: 
Analytical Predictions for Tensing 


Simon Foreman and Leonardo Senatore 

Stanford Institute for Theoretical Physics and Department of Physics, 
Stanford University, Stanford, CA 94306 

Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford University and SLAC, Menlo Park, CA 94025 


Abstract 

We study the prediction of the Effective Field Theory of Large Scale Structures (EFTofLSS) for the 
matter power spectrum at different redshifts. In previous work, we found that the two-loop prediction 
can match the nonlinear power spectrum measured from A-body simulations at redshift zero within 
approximately 2% up to /c ~ 0.6/iMpc“^ after fixing a single free parameter, the so-called “speed of 
sound”. We determine the time evolution of this parameter by matching the EFTofLSS prediction to 
simulation output at different redshifts, and find that it is well-described by a fitting function that 
only includes one additional parameter. After the two free parameters are fixed, the prediction agrees 
with nonlinear data within approximately 2% up to at least k ~ l/iMpc“^ at z > 1, and also within 
approximately 5% up to /c ~ 1.2/i-Mpc“^ at z = 1 and k ~ 2.3/iMpc“^ at z = 3, a major improvement 
with respect to other perturbative techniques. We also develop an accurate way to estimate where the 
EFTofLSS predictions at different loop orders should fail, based on the sizes of the next-order terms that 
are neglected, and find agreement with the actual comparisons to data. Finally, we use our matter power 
spectrum results to perform analytical calculations of lensing potential power spectra corresponding to 
both CMB and galaxy lensing. This opens the door to future direct applications of the EFTofLSS to 
observations of gravitational clustering on cosmic scales. 
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1 Introduction 

The WMAP and Planck satellites, which have measured the properties of the CMB with exquisite 
precision, have enormously improved our understanding of the primordial universe. But they have 
also made any further progress in our knowledge of the initial conditions of the universe terribly 
hard. The next leading source of cosmological information is most probably going to be the Large 
Scale Structures (LSS) of the universe. But in order for them to be able to improve our knowledge 
of the primordial universe, a better understanding of their behavior is required. The Effective 
Field Theory of Large Scale Structures (EFTofLSS) [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15] 
is a research program that aims at developing an analytic understanding of the clustering of dark 
matter and galaxies at large distances, in a perturbative expansion in the smallness of the size of 
the nonlinearities. 
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In LSS, short distance nonlinearities are very large, and therefore not under perturbative 
control. Mode-mode coupling implies that short distance physics does affect long distances, where 
fluctuations are small and should naively be amenable to a perturbative treatment. In order to 
correctly encode the effect of short distance nonlinearities, in the EFTofLSS additional terms are 
present in the equations of motion. The form of these terms is highly constrained by symmetries, 
mainly the equivalence principle that controls the interactions with gravity [1, 2, 8], and by the 
nature of the interactions between different species, such as baryons and dark matter [14]. At 
a given order in perturbation theory, only a finite number of terms is necessary, each one with 
its own numerical prefactor, that we can call “coupling constant” [2, 4, 5]. Therefore, since 
the functional form of each term is fixed, and we are free to adjust a finite number of coupling 
constants, the EFTofLSS is a predictive theory. In fact, it falls in the same class of predictivity of 
general relativity, the Chiral Lagrangian, or dielectric materials. All of these theories are valid in 
a specihc energy range and make predictions, at a given level of accuracy, that depend on some 
unknown coupling constant that needs to be measured in some specihc observations, after which 
all the remaining observations are predicted. In a modern language, all of these theories are called 
Effective Field Theories (EFTs). 

To make the analogy even more concrete, let us focus on the case of dielectric materials. 
When we consider propagation of electromagnetic waves in a dielectric material, we do not solve 
the vacuum Maxwell equations with an Avogadro number of atoms; instead, we solve the dielectric 
Maxwell equations with no sources. The dielectric equations are valid in the range of scales from 
the size of the dielectric material to the atomic scale. Therefore, waves propagating in the medium 
must have a wavelength that falls into this regime. Similarly, the EFTofLSS describes huctuations 
in LSS quantities with wavelengths that are longer than the nonlinear scale. 

In the dielectric material, at a given order in perturbation theory, the effect of the detailed 
structure of the materials can be encapsulated in a few parameters: conductivity, polarizability, 
etc. Similarly, in the EFTofLSS, at a given order in perturbation theory we have to specify a 
few parameters, such as the speed of sound of the species (or more generally the coefficients of 
the various terms in the effective stress tensor) [2, 5, 14], the bias parameters [12], and the short- 
distance velocity contribution for redshift distortions [13] . Back to the case of dielectric materials, 
the actual numerical value of the conductivity can be measured directly in an experiment, or can 
be measured a priori by using a computer simulation. Similarly in the EFTofLSS, the various 
parameters that are necessary to match observations can be measured directly in observations (that 
is without any need of simulations), or, even better, they can be anticipated by measurements in 
simulations. Finally, as the theory of dielectric materials is predictive and can be verihed (and 
it has been verihed) experimentally, similarly, the EFTofLSS is predictive and can be verihed by 
observations. 

So far, the EFTofLSS has been used to predict the equal-time dark-matter power spectrum at 
one- [2] and at two-loop orders [5, 9], the equal-time momentum power spectrum at one loop [9], 
the equal-time dark matter bispectrum at one loop [9] , the slope of the equal-time velocity vorticity 
power spectrum at leading order [5, 6], the equal-time baryon power spectrum at one loop [14], and 
the ehect of baryons relative to dark matter ehectively at two loops [14]. All of these quantities 
have been computed at redshift z = 0. Only one parameter, the so-called ehective speed of sound 
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of dark matter, has been used for the predictions involving just dark matter, and an additional 
one, the so-called effective speed of sound of baryons, has been introduced when dealing with 
baryons. These parameters have been chosen to match numerical observations of power spectra and 
bispectra, and the EFTofLSS has then been able to match the numerical data to percent accuracy 
at /c ~ 0.3hMpc“^ for calculations done at one-loop order, and at /c ~ 0.6hMpc“^ for two-loop 
calculations. In the case of the dark matter speed of sound, this parameter has been also measured 
directly using the dark matter particles as degrees of freedom (this is the analogous procedure 
as predicting the conductivity of a material by solving the quantum mechanical equations of 
atoms), and the results agree with the number obtained by directly matching the long wavelength 
observables such as the power spectrum, within error bars [2], 

These results form a growing body of evidence for the correctness and viability of the EFT 
approach to structure formation. They represent a major improvement for the number of modes 
that are amenable to analytic control. Former analytic techniques, when implemented without 
the addition of extra ad-hoc (that is, physically unjustihed) htting parameters, fail at about 
k ~ 0.1hMpc“^. Since the number of modes scales as the cube of the maximum wavenumber, 
these results suggest that there is a factor of 200 more modes amenable to analytic techniques than 
previously believed. This could have drastic consequences for the amount of physical information 
that LSS surveys can deliver. 

The above statement about the number of modes is made by extrapolating to higher redshifts 
the results obtained at redshift zero, a procedure that is, strictly speaking, unjustified. The 
purpose of this paper is to explore the predictions of the EFTofLSS at all redshifts. Since in the 
universe time-translations are spontaneously broken, the coefficients of the EFTofLSS are time- 
dependent. In fact, since the short modes evolve on the same time scale as the long modes, the 
EFTofLSS is actually local in space, but non-local in time [5, 7]. Therefore, the coefficients of 
the counterterms of the EFTofLSS are integrals over time of some unknown function. Thanks to 
the fact that the linear evolution of the modes is /c-independent, the ignorance of this function of 
time reduces to having a different time-dependent coefficient for each order in perturbation theory 
at which the counterterm is evaluated [5]. In the local in time approximation, these different 
coefficients become equal, and it was shown in [5, 12] that this approximation appeared to be 
preferred by the data. By studying the prediction of the theory as a function of redshift, we will 
explore this time dependence as well. 

Our main avenue for this exploration will be the equal-time matter power spectrum. After 
a brief review of the formalism behind the EFTofLSS, we will describe a method for estimating 
the contributions of various loop orders to the power spectrum prediction, properly accounting 
for the strength of mode couplings at different scales and for important combinations of loop 
diagrams that occur at different orders. We will then embark on a comparison of the EFTofLSS 
prediction to power spectra measured from Wbody simulations at different redshifts. Ending that 
the prediction’s range of validity is consistent with our estimates at all times. Motivated by the 
picture of a scaling universe with a running slope, we present a two-parameter fitting function 
for the time-dependence of the effective speed of sound, and demonstrate that it can be used to 
reproduce the results obtained by fitting the speed of sound separately at each redshift. 

With these results in hand, it is then possible to perform an analytical calculation of the 
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lensing potential power spectrum Cf, and by extension many other observable quantities related to 
gravitational lensing. We present computations corresponding both to lensing of CMB photons and 
photons emitted by galaxies at lower redshifts, and quantify the accuracy of these computations: 
for CMB lensing, the two-loop EFT prediction for Cf has less than 5% theoretical uncertainty 
for i < 1000, while for photons emitted from sources at z = 2 or z = 1, similar accuracy is 
achieved for £ < 600 and 350 respectively. We also estimate the effects of baryonic physics on 
these predictions, finding them to be at percent level in the above multipole ranges. 

2 Review of the EFTofLSS for Dark Matter 

In the EFTofLSS, the dynamics of collisionless dark matter on large scales are described by the 
following Fourier-space equations of motion for the overdensity {6 = Sp/p) and velocity divergence 
(6 = diV^) helds [3] 

oH daS{k, a) -f 6{k, a) = — a{q, k — q)6{k — q, a)6{q, a) , (1) 

Jq 

oH da9{k, a) + ^{k, a) + ^'H‘^{a)£l^{a)5{k, a) = — f f3{q, k — q)9{k — q, a)9{q, a) — iki{dT)p^'‘(k) , 

^ J<5f 

( 2 ) 


where 




nl 


aH?{a) 


a(fci, k2) = 


(fci ^2) ■ ki 


f3{ki,k2) = 


(fci ^ 2 ) ki ■ ^2 


( 3 ) 


and we have used the shorthand notation f^ = f The operator (9r)p;* incorporates nonlinear, 

short-distance physics whose effects on the long-distance modes we would like to include in the 
theory. It is defined by = p~^djT^^, where is the stress tensor, but from now on we 

will use the term “stress tensor” to refer directly to the composite operator We write this 

quantity as an expansion in derivatives and helds, consistent with the equivalence principle, which 
is the relevant symmetry at the scales we wish to describe. 

One may worry that (5r)p;®, which is generated by short modes of the helds that have been 
“integrated out” of the theory could exhibit some degree of non-locality in time, since the 
(nonlinear) short modes will evolve on roughly the same timescales as the (perturbative) long 
modes [3, 7]. However, it was recently argued in [12] that this non-locality might only ahect the 
power spectrum through a correction of order (tshort/^iong) x A-ioop ~ -Ps-ioop, where in fact the 
short modes that contribute most strongly are expected to have tshort/^iong ~ 0(1/10). While this 

^Vorticity is generated in the EFTofLSS [1], but it contributes only at high order [3, 6], and can be neglected 
at the perturbative level at which we work in this paper. 

^In the EFTofLSS, modes are “integrated out” by writing the theory in terms of fields that are smoothed over 
some scale, taking the expectation value over the short modes, and then taking the smoothing scale to zero within 
the long-wavelength theory—see [2] for a precise description of the smoothing procedure and the construction of 
the EFTofLSS. 
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argument is approximate and needs to be checked further, it agrees with the fact that, in [3], it was 
found that nonlinear power spectrum data at ^ = 0 are better described by predictions derived 
from the time-local limit of (9 t)p/ than by a simple parametrization of non-locality. Therefore, 
after reviewing how non-locality was handled in [3] , the body of this paper will focus on the time- 
local case. This will be further justihed by App. E, which shows that non-locality in time does not 
signihcantly alter the EFTofLSS predictions for the power spectrum at any redshift we consider. 

In [3] it was found that for a two-loop calculation of the power spectrum to describe nonlinear 
data with percent-level accuracy, it is sufficient to consider the lowest-order term in the expansion 
of (5r)p,®, neglecting terms of quadratic or higher order in fields, and also neglecting higher 
derivative terms: 


= / dT'Ki(r,r')5*aV(r',a;fi[r,r';a;]) = 


da' 


a''H{a 


—K{a,a')d'6{a',Xf[[a,a'-,x]), (4) 


where we have used Poisson’s equation in the second equality to trade 9^0 for 6, and ki and K are 
unknown kernels with support of order one Hubble time. Using Xf[ as the path of a fluid element, 
defined by 


x^[a^a]x\ = X — 


da" 


a"H{a" 


-v{a",x^[a,a"]x\) , 


( 5 ) 


ensures that the equations of motion are diffeomorphism-invariant (i.e. invariant under generalized 
Galilean transformations). With this form of the Euler equation (2) becomes 


aH da9{k, a) + T-LO^k, a) + «) = ~ / (3{q, k — q)6{k — q, a)9{q, a) 




da' 


:^^^A:(a,a')[^(a',a;fl)]fc , (6) 


where [fixa)]/, = f d^xe One convenient choice for the kernel K{a,a') is 


K{a, a') = 


(27r)c 


di) 




C(C + 5/2)(p-h l)'H(a) /(a) a'n{a') 


,.dDi{a') Di{a'y-^ 


da' Di{a)P~'^ 


( 7 ) 


which parametrizes the non-locality in time by a fixed power p. Here, T*i(a) is the linear growth 
factor, normalized to unity at a = 1: 


Oi(a) ^ . C+(o) ^ |sJ„(a)W’(a) T . (g) 

T>+(1) 2 Jq nya') 

■^In [3], the specific choice C = 2 was used in Eq. (7). This ensured that the tree-level counterterm in the 
matter power spectrum, Ptreei the same time-dependence as the one-loop correction Pi-ioop- This choice has 
a negligible effect on predictions made at a single redshift (such as those presented in [3] aX z = 0), but for the 
purpose of relating predictions at different times, it is important to allow for other choices of See [10] for the 
details of how to perform calculations in the C ^ 2 case. 
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The equations (1) and (6) admit solutions of the form 

OO OO 

6(a.k) = , 

71=1 n=l 

OO OO 

9(0, k) = -H(o)/(o) ^|B,(o)]"9<"l(fe) - H(o)/(o) J]|B,(o)]"+<9<">(fe) , (9) 

n=l n=l 

under the assumption that h2m(a) = where /(a) = c^logDi/c^loga. The error incurred 

using these approximate solutions instead of the exact solutions has been shown numerically to 
be very small in the one-loop power spectrum (e.g. [2]), so we will make use of these solutions in 
the remainder of the paper. In an upcoming paper [16], we will re-examine this approximation in 
detail. 

Plugging the solutions (9) back into the equations of motion and writing 

= f ... f {27rfSB{k-q^ -(10) 

JQi “'<?n 

and analogously writing 5^^\k), and 9^^\k) in terms of Gn, Fn, and Gn, one obtains 

recurrence relations for the F and G kernels, which can be found, e.g., in [10]. From there, any 
correlation function of 5 and 9 can be written down up to a specified order. The basic building block 
of these correlation functions is the linear solution 5^^\k), assumed to have Gaussian statistics 
(in the absence of primordial non-Gaussianity, which we will assume in this work) and power 
spectrum 

{6^^\k)6^^\k')) = (27r)3<5D(fc + fc')Pii(fc) , (11) 

which can be obtained from a public code such as GAME [17]. The presence of (5r)p/ in the 
equations of motion will lead to the appearance of counterterms in expressions for correlation 
functions. From a very practical standpoint, these counterterms are the main distinction between 
the EFTofLSS and other perturbative approaches to structure formation. 


3 The Matter Power Spectrum at Different Redshifts 

In this section, we compare the one- and two-loop EFTofLSS prediction for the matter power 
spectrum to data from A^-body simulations. At a single redshift, the one-loop prediction has 
one free parameter, and so does the two-loop prediction, provided that the two-loop terms are 
renormalized via the procedure we will describe below. At multiple redshifts, it might seem that 
this parameter needs to be ht separately at each redshift, but in fact there is a convenient way to 
parametrize the time-dependence that involves only one additional parameter. Sec. 3.1 contains 
the formulas used and assumptions made in our predictions. Sec. 3.2 shows estimates for how far 
in k the one- and two-loop predictions should be expected to reach. Sec. 3.3 describes a procedure 
for renormalizing the two-loop contribution, and Sec. 3.4 contains the actual comparisons to data. 
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For the nonlinear power spectrnm, we use version 2 (“FrankenEmu”) of the Coyote inter¬ 
polator [18] with WMAP7 cosmological parameters h = 0.704, = 0.272, fib = 0.0455, 

fiA = 1 — fim, ns = 0.967, and ag = 0.81. As a conservative estimate of the error in the output 
of the interpolator, we take the nonlinear spectrum to have uniform 2% errorbars for z < 1 and 
k < IhMpc”^, 3% errorbars for z > 1 and k < IhMpc”^, and 6% errorbars for z > 1 and 
k > 1 hMpc”^ App. D contains a consistency check of the EFT prediction against a different 
cosmology (WMAP3). We hnd results that are comparable to what we present in Sec. 3.4. 

Our primary metric for quantifying the performance of various perturbation theory predictions 
will be the wavenumber at which the prediction begins to deviate from the measured nonlinear 
power spectrum by a certain percentage. In particular, we will choose this percentage to be 
our estimation of the uncertainty on the nonlinear spectrum we use, and dehne k^^^^ to be the 
wavenumber at which the L-loop EFT prediction exceeds this uncertainty ®. For the convenience 
of the reader, the uncertainty on the data is summarized in Table 1. 

3.1 Formulas 

The one-loop EFT formula for the power spectrum is given by 

^EFT-l-loop(A;, z) = [Di{z)fPu{k) + [Di{z)]'^Pi.loop{k) + Ptree(^, z) , (12) 

^The most accurate settings of the Coyote interpolator do not accept h as an independent parameter, but rather 
fix its value based on the CMB distance to last scattering. Therefore, the nonlinear power spectra we compare 
to in this section actually have h = 0.7013 instead of /i = 0.704, but this difference (^0.5% of the value of h) is 
insignificant for the purpose of comparing to EFT predictions. 

®Our errorbars are determined in the following way. From our understanding, the output of the Coyote emulator 
has uncertainties that are quoted as 1% for fc < 1 ft, Mpc“^ and z < 1, and 5% in all the other regimes of interest [18]. 
From Figs. 8-9 of the fourth paper in Refs. [18], it seems to us that the region k < IftMpc”^, z > 1 has a smaller 
error bar than 5%, and we take it to be about 2%. In this paper we are interested in understanding when the 
EFTofLSS fails in the UV to about 1%. So, to the error bars of Coyote, we add a 1% error to represent when the 
theory is within 1% of the numerical data. This error should also account for potential unknown systematics in the 
simulations and in the comparison between simulations and the EFTofLSS. At last, we point out that when the 
EFTofLSS begins to fail in the UV, it does so in a very steep way, so that the UV reach for the EFTofLSS that we 
establish later is expected to be quite insensitive to the precise choice of error bars. Similarly, the statement that 
the UV reach of the EFTofLSS is compatible with one from the estimated theoretical error that we derive later is 
quite unaffected by this. 

®We wish to emphasize that our chosen failure criterion is just one of many possibilities. In this work, we 
are comparing to numerical measurements with a certain uncertainty, and so we find it to be natural to use this 
uncertainty as a benchmark for our comparisons with theory. However, for more detailed applications of the 
EFTofLSS to observational data, or to other simulations with different uncertainties, there can certainly be other 
criteria which could be more suitable. Indeed, if predictions of very high precision are needed at some point in the 
future, the associated failure wavenumber of the EFTofLSS will be smaller than what we quote in this work, since 
the effects of loops and counterterms will be larger relative to the required precision. 



2 ; 

k 

failure threshold 

< 1 

<lh Mpc“^ 

2 % 

> 1 

< 1 h Mpc“^ 

3% 

> 1 

> IhMpc"^ 

6 % 


Table 1: A summary of the uncertainty on the numerical data, and therefore of the criteria we use 
to define the “failure” of a given prediction in this work: when the deviation of the L-loop EFT power 
spectrum prediction from the nonlinear measurements reaches the threshold in the rightmost column, we 
consider the corresponding wavenumber, to be where the prediction fails. 


where while the two-loop formula is 


-PeFT-2-1oop(^, 2 ;) — -PEFT-l-loop(fc5 + [Dl{z)]^P2-\oop{k) — 2{27l)c 

+ (2;r)cI(i)(^)l-Di(z)l'‘-Pi'.toh(*^) + ( 2 ! r )2 (^1 + |c2(i) 


^.p,(z)Ppn(k) 

^NL 

{z)riD,{z)fP 

^IMT 


Pu(k) . 


(13) 


Expressions for Pi-ioop{k) and P 2 -ioop(fc) are given in [5], and also in [3]. Expressions for -Pi 4 oop(^) 
are given in [3] for the C = 2 case (see Eq. (7)); the extension to ^ 7 ^ 2 is straightforward using the 
recurrence relations from [10]. The coefficient of the {k/k-!<si^)^Pii term is different from that given 
in [3, 9], owing to the inclusion of a tree-level counterterm that has previously been omitted from 
two-loop computations; the derivation of this term (along with the other tree-level counterterms 
shown above) can be found in App. A. 

The time-dependence of comes from the fact that the diagrams in are built in one 

of two ways: by multiplying Pj^^fl/Pn, which scales like c^(i)(- 2 ), by a diagram contained in Pi_ioop, 

which scales like [Di{z)Y] or by inserting P^^l/Pn into a one-loop diagram. This second type of 
diagram may have time-dependence that could differ slightly from c^^^.^-^{z)[Di{zyf, determined by 
the detailed form of the integrals over Green’s functions in time. However, just as the Green’s 
function integrals in Pi_ioop collapse to ~ [Di^zff to very good accuracy in AGDM, we expect 
that the integrals in will also collapse to c^^f^^-^{z)[Di{zyf at a similar level of accuracy, and 

therefore we will use this expression in our calculations. Furthermore, the Pifioop (fc//cNL)'^Pii 
terms depend slightly on the value of C, we choose, so a particular value must be chosen before 
the prediction is computed; we use ( = 3, which is very roughly motivated by scaling universe 
arguments '. Finally, we use the form of Pi^oop corresponding to the time-local case (the p ^ 00 
limit of Eq. (7)). In App. E, we show that relaxing the assumption of locality in time does not 
signihcantly affect the comparison to data. 

^In principle, the time-dependence of and therefore the value of C, can be measured from fitting to sim¬ 
ulations. However, as just mentioned, it is useful for the calculation to make an initial assumption about the 
time dependence. The fitting to the data that we will explain and perform later will be such that the residual 
dependence on this initial choice is very small. 
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Figure 1: Left: effective slope nes{k) = d\ogPrm{k)/dlogk, where Pnw is the linear matter power 
spectrum without BAO wiggles. Pnw(^) is given by Eg. (If)- Center: running of this slope, dnes/dlogk. 
Right: running of the running, d^ne^/dlogk'^. The large running of the slope implies that a simple 
analogy with the case of a pure scaling universe will he insufficient to approximate the behavior of various 
terms in the power spectrum prediction. 


For numerical computations, we use CAMB [17] to generate a linear power spectrum, and 
a modification of the Copter library [19] that makes use of the IR-safe integrands from [5] and 
Monte Carlo integration routines from the CUBA library [20] to compute the loop integrals. 


3.2 Estimates and expectations 

An essential feature of the EFTofLSS is that it is possible to estimate at what scale a given 
prediction should fail, based on the size of the next-order terms that are neglected. In [3, 10], 
these estimates were performed by treating the predictions as taking place in a scaling universe 
{Pii{k) oc kP) with slope n determined from the slope of the ACDM linear power spectrum around 
the scales where the two-loop z = 0 prediction should be valid {k < O.OhMpc”^). However, if 
one is interested in a wide range of scales and higher order contributions, where errors are given 
the possibility of piling up, the scaling-universe analogy can no longer be simply applied. In fact, 
in Figure 1, we plot the effective slope nes{k) of the linear power spectrum without BAO wiggles, 
as found from the fitting formula in [21] as well as the running of this slope, dries/d\ogk, and 
the running of the running, d'^Ues/d log k"^. While d^Ues/d log k"^ is quite small, the running is not 

®For reference, this formula is given by 


P„w(fc) cx 


log(2e-|- 1.8g) 


where 


_log(2e -f 1.8q) + q‘^{U.2 + 731/(1 -f 62.50). 

Teff = “t“ 

which itself depends on 


q = 


llbsVr-i 


/iMpc"^ J \ 2.7K J 


eff ’ 


1 -f (0.43fcs)^ 


4 / ’ 


(14) 

(15) 


44.5 log (9.83/Mm/ifi) , Bb , 9 , ( 

- -hi- .A Mpc , ar = 1 - 0.3281og(331Bmh07^ -t 0.381og(22.3BmhO • (16) 

[l + 10(Bbh2)F4]^/'' 


(The prefactors of Pnw are not needed to calculate neff(fc) = dlogPnw(fc)/dlogfc.) 
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small, at least for k < 0.6/iMpc“^, which means that the universe is not well approximated by 
a scaling one. Furthermore, and maybe more importantly, there can be important cancellations 
among diagrams that are not captured by the simple scaling argument, as we will now describe. 
This pushes us to develop a method to estimate the various contributions in a way that mirrors 
more closely the perturbative structure of the theory. 


3.2.1 One loop 


Let us begin by formulating such an estimate for Pi_ioop- To do so, we recall that Pi-ioop primarily 
captures the effects of density fluctuations from modes with q < k and displacements from modes 
with g > A; on a density mode of wavenumber k [8, 9]. These effects are respectively encapsulated 
in the parameters es< and es>, dehned by ® 

pk 2 3 

e5<ik)Pii{k) = dgpi_ioop(fc,g)|g«fc Pii(A;)-— / dgg^Pii(g) , (19) 

Jo Jo 

1 / ^22k‘^\ 

-es>{k)Pii{k) = dqpi.ioop{k,q)\g^k-Pii{k)j^ 

We therefore use 

d‘ro“’W“(«i<-£.>)niW (21) 

as our estimate for Pi_ioopik) = Pi-ioop(A)) + Pl;^g{k). (This estimate assumes not more than an 
order one cancellation between Pi_ioop and at all redshifts.) In order to use this estimate 
without knowledge of the exact value of the parameter that we will hnd by comparing to 
data, as well as of order one numerical factors resulting from the loop integrals, we multiply 
and divide the right-hand side of (21) by 2 and use the results as a range for the estimate of 
^i-ioop- indicator of where linear theory should fail (i.e. as dehned at the beginning 

of Sec. 3), we hnd where our estimate for the ratio Pi_*io*op^(^)/Pii(A:) equals the uncertainty on 
the nonlinear spectrum. The result is that we expect the failure to occur somewhere in the 
range 0.04hMpc“^ < kj^^ < 0.07hMpc“^ at z = 0, while the actual failure occurs at = 
O.OShMpc”^. Fig. 2 displays the estimated and exact failures at various redshifts from 0 to 4. 

®More in detail, the parameters upon which the Eulerian-space EFT Taylor expands are [8, 9] 


€S< 



(27r)3 


Pii(P) , 


,2 r Pn(fc') 

A (2^)3 


P Pk’ Pii(P) 
A ( 27’‘)3 fc'2 


(17) 


They represent the effect of long-wavelength overdensities, the effect of short-wavelength displacements, and the 
effect of long-wavelength displacements. Here long and short are relative to a specific wavenumber k that we are 
interested in. In the case of a system with two species, such as dark matter plus baryons [14], we also have 


rel 

^s< 


{k) 


Pk' Pii(fc') 

(27r)3 fc'2 > 


rel 

^s> 


{k) 


2 /•“ Pk' Pn(fc') 
A (27r)3 k'^ 


(18) 


where Pn is the power spectrum of d{Sc — Sp /dloga, and are the overdensities of the dark matter and 

baryons respectively. After IR-resummation [9, 14], the parameters and are not expanded upon. 
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Linear theory 


z 

estimated failure 
[h Mpc“^] 

0 

0.04-0.07 

1 

0.06 - 0.09 

2 

0.08 - 0.17 

3 

0.11 - 0.23 

4 

0.13 - 0.31 


Figure 2: Left: Linear theory prediction for the power spectrum, normalized to the nonlinear spectrum 
from the Coyote emulator at different redshifts. The dashed lines show the estimated error on the nonlinear 
data, while the dotted line at the value Pii/Pcoyote = 1 is provided as a visual aid. Right: Estimates for 
when linear theory should fail, based on when the ratio P[[°lfp /Pii exceeds the error on the nonlinear 
data. We quote a range of wavenumbers based on dividing and multiplying the estimate from Eq. (21) 
by a factor of 2. By comparing these estimates with the plots on the left, we find that linear theory fails 
roughly where it should. 


For a more precise estimate, instead of multiplying and dividing Eq. (21) by 2, we could add 
the Ptreg(/c) term, using the exact value of obtained from one-loop hts to the data. Doing this 

gives = 0.06hMpc“^ at z = 0, and also gives estimates at other redshifts that are similarly 
close to the actual failure of linear theory. We will follow this kind of strategy when estimating 
P^^ioop cind sums of all terms entering at two- and three-loop order, respectively), in 

order to obtain the most precise picture of when the one- and two-loop EFT predictions should 

fail. 

3.2.2 Two loops 

To estimate the size of P 24 °*op^ we make use of the fact that P 2 -I 00 P will be dominated by diagrams 
containing two independent trivial angular integrations, which each contribute a factor of (27r). 
The sum of these diagrams can be approximated, up to combinatorics, as a product of the standard 
one-loop diagrams, Pi_ioop, and a P 13 diagram with the Pii(fc) factor removed (see Fig. 3).^° 

^*^When the diagrams are written using only cubic vertices instead of symmetrized n-point vertices (e.g. [16]), 
there will be also diagrams with two independent angular integrations but whose loops are nested. However, these 
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Figure 3: Schematic depiction of the dominant two-loop diagrams in the power spectrum: [P 22 + -P 13 ] x 

p(no-IR) ^ p_^_^ 

However, an exact calculation of Pi 3 {k) will include a contribution from modes with momentum 
q k, which will cancel with the corresponding contribution from other diagrams (at one-loop 
order the cancellation is provided by P 22 {k)). We should therefore include only the part of P 13 , 
that we call that is not going to cancel against other diagrams in the hnal resnlt. This 

contribntion cannot be calculated precisely, but a rough estimate can be obtained by hnding 
where the P 13 and P 22 integrals exhibit strong cancellations when each integral in wavenumber q 
is cut off at some gmax, and then taking P{^ as the contribntion to P 13 from q < gmax- We hnd 
that gmax is somewhere between /c/5 and 2/c/5 for k < 0.4/iMpc“^, and this leads to = 

Pu{k) - P,T\k) — cx'Pi^iJf) with 1/10 Oi' "S 1/3. Also, the quantity that appears in T-loop 
contribntions will not be Pi_ioop or a'Pis on their own, bnt rather the combinations Pi_ioop + -Ptrei 
or q;'Pi 3 -|- Ptree- For convenience, we can use the fact that P 13 ~ 3Pi_ioop in the range we care 
about to replace a'Pis by aPi.ioop in our estimates, so that 3/10 < a < 1. 

In snmmary, then, we estimate 

“<3 [u,o„p(i) + utA)] . (22) 

with ranging between 1/2 and 2 to acconnt for nnknown combinatorics and other order one 
nnmerical coefficients. This expression highlights the possibility that there may be rather large 
cancellations between loops and connterterms at two-loop (and higher) order, because loops be¬ 
yond the first contribnte factors of Q;Pi_ioop + Ptree instead of Pi-ioop + Ptree (which was assumed, 
and fonnd, not to exhibit a cancellation beyond order one). 

We can explicitly see such a cancellation occurring if we try to match our estimate for P 24 °o^p 
with the exact calculation. With f3 hxed to 1, onr estimate matches the shape of the exact P 24 °o^p 

diagrams will contain more integrations over Green’s functions in time. Recall that in the Einstein de-Sitter case, 
these Green’s functions are roughly the scale factor a, so m integrals over Green’s functions will give an overall 
prefactor of ~ 1/m!. This factor will suppress diagrams with nested loops relative to those with non-nested loops, 
so we only consider the latter type of diagram in our estimates. 
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z = 0 



■ 0 . 02 ^.- 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


k [h Mpc ' ] 

Figure 4: Exact and approximate versions of P 2 °imp /-^n at z = t), the latter plotted using Eq. (22) with 
a ~ 3/10 and (3 = 1. The bottom panel is a zoomed-in version of the top panel. Overall, the shapes of 
both curves match very well up to k ^ 0.6/iMpc“^ . 

if a ~ 3/10 (see Fig. 4) Examination of after has been matched to data reveals 
that ~ —(l/3)Pi_ioop at 2 ; = 0 and k < 0.4hMpc“^. Therefore, at these values of and 

a, the factor (aPi_ioop + Aree) (22) heavily suppresses the magnitude of P 24 °*op^ in this range. 
This suppression was hrst noticed in [3] as a cancellation between the different terms entering 
the power spectrum prediction at two loops, but it is more naturally understood as arising from 
the relationship between loops and counterterms when they are present as subdiagrams within 
higher-loop diagrams. 

When attempting to estimate P 24 °*/]|^ at higher redshifts, we hnd that some choices for a imply 
that the expression (22) will have zero crossings that are not expected to be present in the full 
calculation, or, if present, clearly cannot be exactly reproduced by our estimates (see Fig. 5 for 
an example aX z = 1). This is particularly problematic because we would like to use the prefactor 
(3 to incorporate the uncertainty in the overall amplitudes of the estimates, but points close to a 
zero crossing will not be rescaled in the desired way if (3 is rescaled. Therefore, at redshifts where 
these zero crossings appear in Eq. (22), we will modify the estimate to preserve its fc-dependence 
at low k when a ~ 3/10, but reflect it across the F-axis so that the estimate no longer changes 
sign (see Fig. 5). Shifting the value of a is the minimal way to accomplish this goal. At the 
redshifts where there are zero crossings, namely z = 1,2, 3, we hnd that the requisite shifts are 
Aa ~ —1/6, —1/3, and —1/2 respectively. 

To obtain a rough value for the wavenumber at which Peft-i-Ioop(F) will fail, we will use the 

similar estimate for the contribution of the two-loop finite term is presented in App. B.l, while comparisons 
between the estimated and exact versions of P 2 -I 00 P different redshifts can be found in App. B.2. 
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p(total) j p 
: —^2-\oop/^n 

estimate with a = 3/10 / : 

p(total) j p 
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estimate with a = -1/6 + 3/10 / 
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/fc[/!Mpc-'] 


Figure 5: Approximate versions of P 2 -°oop /at z = 1, plotted using Eq. (22) with a ~ 3/10 (dashed 
curve) and —1/6 + 3/10 (dotted cruve). Both curves use /3 = 1. The dashed curve has zero crossings that 
are not expected from the exact calculation, while the dotted curve has the same slope as the dashed one 
at low k, but reflected across the k-axis, eliminating any zero crossings. This shift in the value of a is 
the minimal way to modify our estimate to remove the zero crossings; we also implement similar shifts 
at z = 2 and 3. 


scale at which the estimate for the ratio i^ 24 oop^ (A;)/Peft-i- 1 oop(^) equals the error on the nonlinear 

data. In our estimate for P 24 °lp^ {k), we allow a and fl to vary within the ranges specihed above: 
3/10 < a < 1 (appropriately shifted at redshifts where there is a zero crossing in the estimate) 
and 1/2 < /3 < 2. We also allow for to deviate by 20% from its best-£t value (taken from our 

hts to the Coyote emulator in Sec. 3.4). At z; = 0, this results in an expectation that should 
fall somewhere within the range 0.14hMpc“^ < k^^^^ < 0.58 hMpc“^, encompassing the actual 
value of = 0.36 hMpc“^. Table 2 shows estimates for for redshifts from 0 to 4, along with 
exact values obtained from the exact hts in Sec. 3.4. 

Before moving on, we remark that P13 not only has an IR contribution that cancels with P22, 
but it also has a UV contribution (from modes with q > /cnl) that is not under perturbative 
control and degenerate with a contribution from a counterterm. Estimating this UV contribution 
at 2 ; = 0 for k ~ 0.4hMpc“^ by integrating Pl 3 from q = 0.7hMpc ^ to 00 , we find that it 
can be as much as 50% of aPi-ioop, indicating that an unrenormalized calculation of a one-loop 
subdiagram can contain a large uncontrolled contribution. In this light, it is not surprising that 
the P/rei term turns out to be a sizable fraction of the size of Pi_ioop after is fit to data, since 

Pt//i must correct for UV contributions to Pi_ioop that should not be included. 

3.2.3 Three loops 

For we follow the same strategy as for P^^ioop, multiplying the lower-loop contribution by 

a factor representing the additional loop: 

pho“;’w«/3 


Q^-Pl-loop(^) ^tree {3. 

Pii{k) 




(23) 
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One-loop EFT 

Two-loop EFT 

2 ; 

estimated failure 
[hMpc“^] 

actual failure 
[h.Mpc“^] 

estimated failure 
[h.Mpc“^] 

actual failure 
[hMpc“^] 

0 

0.14 - 0.58 

0.36 

0.26 - 0.78 

0.65 

1 

0.20 - 0.60 

0.45 

0.42 - 0.85 

1.2 

2 

0.27-0.85 

0.68 

0.65 - 1.50 

1.6 

3 

0.5 - 1.4 

1.15 

1.2 - 3.0 

2.3 

4 

0.5 - 2.3 

1.65 

1.3 - 5.0 

3.3 


Table 2: Estimates for where the one- and two-loop EFT predictions should fail, based on where the 
ratios P 2 -°ioop I^^PT-i-ioop and / Peft-2-Ioop respectively exceed the error on the nonlinear data we 

use. We also state the scales at which the exact one- and two-loop EFT calculations fail, as seen in the 
comparison to data in Sec. S.f. The details of the estimates can be found in the main text. The exact 
two-loop calculation sometimes reaches slightly farther than the estimates suggest, but this is not cause 
for alarm given the 0(1) uncertainty in the estimates. Note also that there is some sizable uncertainty 
in the determination of the actual failure points, stemming from the uncertainty in and the range of 
reasonable choices for kreni see the main text for further discussion. 

We use the exact calculation for and also use the same ranges for a, (3, and 

as before. Similarly, the estimate for where Peft- 2 - 1 oop should fail, denoted by is taken to 
be where ^’ 34 °opV-PEFT- 2 -ioop exceeds the error on the nonlinear data. At 2 ; = 0, Eq. (23) gives 

0.26hMpc“^ < < 0.78hMpc“\ which includes the exact value of = 0.65 h Mpc“^. 

Values for other redshifts are shown in Table 2 . At some redshifts these estimates fall slightly 
short of where the exact calculation fails, but this is not cause for alarm, given the fact that 
these estimates are only intended to be accurate up to 0(1) factors. The fact that the estimated 
and exact reach of Peft- 2 - 1 oop show the same trend indicates that the suppression due to the 
(aPi_ioop + -Ptoee) factor that we discussed at two loops indeed persists at higher loops, since 
otherwise the prediction would fail earlier. 

Note that there is some sizable uncertainty in the determination of the actual failure points 
listed in Table 2, stemming from the uncertainty in and the range of reasonable choices for 
/cren; this is shown by the blue bands in Fig. 8 . Nt z = 0, this uncertainty could result in an actual 
failure point of around k ~ 0.2 h Mpc“^ for Peft- 2 - 1 oop, because of the sensitivity of the low-redshift 
calculation to accidental cancellations between different terms, k ~ 0 . 2 hMpc“^ is slightly below 
the lower bound of our z = 0 estimate in Tab. 2, but this is not cause for concern since the 
estimates themselves have 0(1) uncertainties. Meanwhile, at higher redshifts, the uncertainty in 
the failure point of the prediction is much less severe, and the range of possible failure points 
indicated by Fig. 8 lies comfortably within our estimates. 
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Figure 6: Ratio of Peft-2-Ioop/Peft-i-Ioop at various redshifts, with set to what one would get 

from a fit of Peft-i-Ioop to nonlinear data, and c^^ 2 ) i^) ^^t to a value that keeps the ratio flat for as large 
a range as possible at low k. (At z >1, this translates into having the ratio as close to unity as possible, 
while at z = 0, the flat region is offset from unity by 1% —see the main text for further discussion.) 
The highest scale where the ratio exceeds a 1% deviation from unity (shown in dashed lines) determines 

hren{,z). 


3.3 Determining c^^ 2 ) 

In moving from the one-loop to the two-loop prediction for the power spectrum, we must include 
several new terms, the most important being the term containing the two-loop diagrams them¬ 
selves, P 2 -ioop- This term will have a contribution that scales like {k/Pn, but this contribution 
arises from UV modes that are not under perturbative control, and whose impact on the power 
spectrum is degenerate with the choice of Therefore, this contribution must be cancelled by 
adding a counterterm ~ c^^ 2 )(^/^nl)^-Pii and adjusting c^^ 2 ) appropriately. We emphasize that 
adjusting the coefficient of the {k/k^EYPu term in this way does not change the value of the 
physical parameter rather, the only role of c^( 2 ) is to cancel the unphysical part of P 2 -I 00 P 
that is erroneously included in the calculation. 

In a scaling universe, one could analytically calculate the form of c^^ 2 ) Taylor-expanding the 
integrand of P 2 -I 00 P for large internal momenta (or alternately for small k): integrating the result 
will give a constant (dependent on the method used to regulate the integral, such as a hard UV 
cutoff A) times kflPn, and c^( 2 ) can then be set to to cancel off this term. In Appendix C, we argue 
that in ACDM, it is dangerous to use this procedure if one does not have access to nonlinear” 
power spectrum data that are very precise at low k. The data we currently use do not have the 
desired precision, and therefore we must use a different strategy for determining our case. 

In [3] and [9], the prediction was only considered at a single redshift, and c^^ 2 ) ''^us fixed by 
simply requiring that the one-loop and two-loop EFT predictions matched either over a certain 
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range of k or at one specific renormalization scale k-^en- We will adopt the latter approach here, 
bnt to do so, we mnst hrst dehne a prescription for choosing fcren at each redshift. In general, fcren 
should be chosen in the range where -P 2 aoop^ small, a fact that pushes fcren —^ 0, but also where 
the k‘^Pii{k) contribution of P 2 -I 00 P is niuch larger than the error in the numerical data, pushing 
/Cren “^ cxD. From a practical standpoint, this will occur where the one-loop and two-loop EFT 
predictions begin to depart from each other; more specifically, in this work we set fcren to be the 
scale at which this departure exceeds 1 %, when is fixed by a fit of Peft-i-1oop to the nonlinear 
data and ^ value that keeps the ratio -Peft-2-1oop/-Peft-i-1oop hat for as large a range 

as possible at low k. Several examples of the ratio -Peft-2-1oop/-Peft-i-1oop are shown in Fig. 6 , from 
which /cren at each redshift can simply be read off. 

Note that at 2 ; > 1, when we apply the procedure described above, Teft- 2 - 1 oop/-Peft-i- 1 oop is very 
close to unity at low A;, while ai z = 0 the ratio is systematically offset from unity before deviating 
at greater than 1%. We allow for this deviation at low z to account for possible numerical error in 
our loop integrals or IR-resummation. Alternatively, if we are very strict with the low-A; flatness 
criterion, we instead find that the ratio deviates at 1% starting around A: ~ 0.1 hMpc“^, leading to 
kren = 0.1 h Mpc“^ instead of 0.3 hMpc“^ at z = 0, where however the measurement of would 
be challenging as its effect is dangerously closer to the numerical errors. We defer a more thorough 
investigation of these two options to future work. In this paper, we use k^en = 0.3hMpc“^ at 
z = 0, but incorporate the uncertainty in k^en into our final results, as described in the next 
section. 

3.4 Comparison with simulations 

We are now ready to compare the two-loop EFT prediction for the matter power spectrum, 
Eq. (13), to results from simulations at different redshifts. (In fact, we will use the IR-resummed 
version of (13), the details of which can be found in [9].) There is only one free parameter, 
Cs(i)(^), to bo £t at each redshift, since the procedure described in Sec. 3.3 fixes c^( 2 )(^) once 
has been determined. At each redshift, we determine a best-£t value of by the 
value for which the ratio TEFT- 2 -ioop/-Pcoyote is constant over as large a range of wavenumbers as 
possible. As a conservative estimate of the uncertainty in this value, we vary in either direction 
until the ratio TEFT- 2 -ioop/-Pcoyote deviates from its constant value by 3% at the upper end of the 
constant region. We apply this procedure at each redshift used to build the Coyote emulator, 
2 ; G {0,1/9,1/4, 3/7, 2/3,1, 3/2, 2, 5/2, 3,4}. Furthermore, we incorporate our uncertainty in the 
optimal value of k^-e^ by repeating this procedure at various A;ren values that are smaller than the 
value determined in Sec. 3.3 

The resulting values are shown in Fig. 7, divided by [Zli( 2 :)]^ to decrease the range of 

^^For 0 > 2/3, we consider /cren values as much as 30% below the value we choose to calculate the main results 
for the two-loop prediction (we do not allow Aren to range above the value from Sec. 3.3, since it is likely that this 
would lead to significant contamination of our matching procedure by finite two-loop terms). For z < 3/7, we allow 
Aren to range down to where the flatness criterion on the ratio FEFT- 2 -ioop/Rc;oyote is most strictly satisfied—for 
example, at z = 0 we consider the range 0.1 AMpc“^ < Aren < 0.3AMpc“^. Of course, a lower Aren will lead to 
higher contamination of the measurement of by numerical error. 
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Figure 7: Top: The values of with associated errorbars, found by matching the two-loop EFT 
prediction for the power spectrum to the outputs of the Coyote emulator at different redshifts, in units 
of (27r)“^(A:NL/)^. The specific matching procedure is described in the main text. The red curve 
shows the fitting function (25) with fi = —1.44 (the best-fit value) and higher coefficients set to zero. The 
good match between the red curve and the points demonstrates that Eg. (25) is successful in describing 
the time-dependence of c‘^^^^{z) using only two parameters, c^(i)(0) and (3. Bottom: The corresponding 
values o/c^^ 2)7 again in units of (27r)“^(A:NL/ found from the procedure from Sec. 3.3. 


the plot. The z = 0 value, 





1 


/ ^NL A ^ 
V hMpc“^ / 


(24) 


is consistent with that found in [3, 9] for a similar cosmology, a reassuring sign that minor variations 
in the background cosmology translate into miniscule variations in the value of In a scaling 

universe, c^(^i^{z) would scale with time like a power of Di{z) [ 6 , 11]. As shown in Fig. 1 , the 
effective tilt nes{k) of the linear power spectrum in ACDM deviates signihcantly from the scaling 
form. Nevertheless, a useful starting point for parametrizing the time-dependence of is to 

suppose that it is determined by the effective tilt at /c = /cren, along with corrections arising from 
the running of 


= c2(^)(0)[F)i(z)]3+a^w , N{z) = nes{k)+(3 


dnesjk) (Pnesjk) 

d\og{k) d\og{ky 


+ ■■■ 


, (25) 


k—kren (■2) 


where (3 and 7 are unknown coefficients, and the ellipsis denotes higher logarithmic derivatives 
of ne?i{k). We use k = fcren because it is the scale at which the effects encapsulated by will 

^^However, larger variations in are noticeable when it is measured from cosmologies with different values of 
CTg: in [10], we found c^(i)(0) ~ 2.45 x (27r)“^(fcNL//iMpc“^ in a cosmology with cts = 0.9, while in Appendix D 
we find ~ 1.31 x (27r)“^(fcNL//iMpc“^ in a cosmology with erg = 0.74. Taken together with the 

measured in this section (from a cosmology with erg = 0.81), these points suggest very roughly that c^(i)(0) oc 
although further investigation will be required to verify this behavior. We leave this to future work. 
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affect the power spectrum quite strongly. Apart for the hrst derivative term in (25), the higher 
derivative terms are very small, so it should suffice to consider only the hrst correction to Uefr due 
to running. In fact, Ueff plus the hrst correction is able to match the htted values very 

well, with a best-ht value for f3 of —1.44 ± 0.42 (htting f3 and 7 simultaneously does not alter the 
results). The best-ht curve is shown in red in Fig. 7 

Strictly speaking, the functional form in (25) is in disagreement with the power-law time- 
dependence, c^(i)( 2 ^) C)C with C ~ 3) been repeatedly assumed during the paper 

(as for example in the assumption that the time-dependence of the terms proportional to c^(i)(^) 
is given by the EdS approximation, as well as in the derivation of the expressions for the higher 
order counterterms). It is possible to relax these assumptions by exactly solving for the time- 
dependence of the diagrams, by using the full Green’s functions as for example done in [ 2 ]. We 
have performed this calculation (which will be presented in an upcoming paper [16]), and found 
that the results are virtually unchanged. 

Comparisons of various theoretical predictions to nonlinear spectra at different redshifts are 
shown in Fig. 8 . The solid blue curves show the two-loop EFT prediction with c^(^^( 2 ;) measured 
separately at each redshift, while the purple dot-dashed curves instead use Eq. (25) with (3 = —1.44 
and setting all higher coefficients to zero. For reference, we also show the one-loop EFT and two- 
loop SPT predictions Fig. 9 shows all two-loop EFT curves on the same plot. As previously 
found in [9], the resummed two-loop prediction matches the nonlinear power spectrum at z = 0 
within 2% (the estimated uncertainty on the nonlinear data) for k < 0.6hMpc“^ after fitting a 
single parameter, that can be determined from data at a single redshift. By fitting at 
each redshift, or alternatively making use of Eq. (25) with only one extra free parameter, it is 
possible to match the nonlinear spectrum at all redshifts up to 2 ; = 4 and up to scales consistent 
with our theoretical expectations of how far the prediction should reach (see Sec. 3.2). 

One can quantify the improvement of the two-loop EFT prediction with respect to two-loop 
SPT by the increase in the number of modes reliably described by each prediction, as determined 
by the scale /cfaii at which each prediction deviates from the nonlinear data by an amount equal to 
the uncertainty on the data. The number of modes scales roughly like so the improvement can 

be estimated by is also possible to perform a more detailed mode-counting 

calculation, based on integrating the number of modes contained in thin shells between z = 0 
and the redshift of interest The results of both calculations are shown in Table 3 for several 
redshifts. The gain in the number of modes is necessarily smaller at higher redshift, since the 

^^Note that a very good approximation for the k^eniz) values we have used in this section is k^eniz) = (0.3 -I- 
0.03z^) /iMpc”^. Using this formula, Eq. (25) can be computed as a smooth, analytical function, which may be 
useful for readers interested in using the EFTofLSS power spectrum prediction in other applications. 

^®We compare with SPT and not with other techniques, such as RPT or LPT, because these and all other 
techniques who differ from SPT only by a resummation of the IR-modes are such that, if they are correctly 
formulated, they must give exactly the same UV reach for the dark matter power spectrum as in SPT [9]. 

^®The details of this calculation are as follows. Starting from the comoving distance between an observer and a 
point at scale factor a. 
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Figure 8: Various theory curves, normalized to nonlinear power spectra at five different redshifts. The 
black dotted lines show the estimated uncertainty on the nonlinear spectra, while the blue band shows the 
uncertainty on the two-loop EFT prediction based on the uncertainty on at that redshift. Whether 
c^(i)( 2 ;) is determined separately at each redshift (blue solid curves) or parametrized using Eq. (25) (purple 
dot-dashed curves), the two-loop EFT predictions agree with the data to much higher k than two-loop SET, 
in a manner consistent with the estimates of the theoretical error from Sec. 3.2. 


the volume of a single spatial shell between scale factors a and a + da is given by 

Fsheii(a) = Wx^{a)dx = 47rx^(a)^da = . (27) 

da a^H (a) 

The number of modes within the volume between redshifts zero and z is found by integrating the number of modes 
within each shell, estimated by the volume of the shell divided by the “fundamental cell” corresponding to fcfaii at 
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Figure 9: The (resummed) two-loop EFT predictions for the matter power spectrum at five different 
redshifts, displayed on the same set of axes and normalized to nonlinear spectra at the same redshifts. 



/l(EFT) /DSPT)Y 

/%ail J 

ivS/NS 

At(EFT) 

modes 

1/9 

(0.66/0.08)^ = 562 

610 

1.5 X 10^ 

3/7 

(0.88/0.15)=^ = 202 

246 

1.1 X 10^ 

1 

(1.2/0.35)=^ = 40 

34 

2.4 X 10® 

2 

(1.6/0.45)=^ = 45 

44 

2.2 X 10^ 

3 

(2.3/0.5)=^ = 97 

59 

8.2 X 10^ 

4 

(3.3/0.8)=^ = 70 

71 

2.4 X 10^° 


Table 3: The ratio of the number of modes reliably described by the two-loop EFT and two-loop SET 
predictions for the matter power spectrum, calculated in two ways: a simple estimate using the fact that 
-^modes oc and a more detailed estimate based on integrating A^modes over thin redshift shells. We also 
show the absolute number of modes described by the two-loop EFT prediction. Depending on the redshift 
ranges of interest, the EFT allows a gain of anywhere between ^30 and ^600 times as many modes as 
SET. 


power spectrum around the quasi-linear scale approaches the slope of u = —3 where loops become 
more important than counterterms. Nevertheless, the EFT allows one to describe a factor of ~44 
more modes out to z = 2, or ~ 70 at z = 4, which can in principle lead to a huge gain in the 
information we can extract from ongoing and upcoming galaxy surveys. 

Finally, let us briefly comment on our procedure for determining at each redshift. An 
alternative procedure would be performing a least-y^ £t over a certain range of k, but there are 
a few obstacles to implementing this procedure in practice. One is that the range over which one 
should £t at each redshift is quite uncertain—the estimates in Sec. 3.2 provide a rough guideline, 
but the resulting best-£t value for will strongly depend on the extent of the chosen range. 

the shell’s redshift: 


N, 


modes 




Fsheii(a) _ A da Wx^jd) 

j(z=o) [27r/fcfaii(a)]^ Ji/{i+z) a'^H{a) [27r/A;faii(a)]^ 


( 28 ) 


22 


















Another obstacle is that there may be systematic errors in the output of the Coyote emulator, 
which would bias our hts even if these errors were within the tolerances quoted in [18]. For 
example, we expect that the “bump” seen at 2 ; > 3 in Fig. 8 is a property of the emulator rather 
than the EFT predictions, since it is also present in T’sPT- 2 -ioop/-Pcoyote- A least-y^ fit would not be 
able to incorporate this piece of information without an extra ad hoc prescription, such as inflating 
the errorbars on the nonlinear spectrum in the neighborhood of the bump. The method we use 
allows us to account for the presence of the bump, at the possible expense of overtuning the fit, 
but we have chosen conservative errors (which incorporate our uncertainty in the renormalization 
scale fcren) cind consistent fc-ranges on our determination of to allow for this possibility. 

Better methods can in principle be designed to measure directly from the power spectrum, 
depending on the quality of available numerical data. Alternatively, can be measured directly 
in simulations using the dark matter particles as degrees of freedom, as done in [2]. 


4 Analytical Calculation of Lensing Potential 

Now that we have the prediction of the EFTofLSS at all redshifts, we are able to calculate various 
observables related to gravitational lensing in a purely analytical way, without relying on sim¬ 
ulations In particular, we will focus on the lensing potential ip, defined so that its gradient 
along the line of sight gives the total angle through which a photon is deflected due to the grav¬ 
itational potential along the photon’s trajectory. If one can reliably compute the statistics of ip, 
one can then straightforwardly compute various observables (such as shear correlation functions) 
that are relevant to weak lensing surveys (see, e.g., [22] for a review). Alternatively, a calcula¬ 
tion of the angular power spectrum of p) can be directly compared to reconstructions from CMB 
measurements. 

To calculate the statistics of pj, we first need the power spectrum of the gravitational poten¬ 
tial (p, defined by 

27r2 

(0(fc, a)p>{kp a)) = —Snik + k')P^{k, a) ; (29) 

this can be related to the matter power spectrum via Poisson’s equation, yielding 

=- —■ (30) 

After expanding the lensing potential in spherical harmonics, tp{n) = ’^i^pJemYemin), we can 
define its angular power spectrum by 

{P^tmpJpm') = . (31) 

In the Limber approximation and assuming a flat universe, Cf is given by a line-of-sight integral 

already stressed, the parameter can in principle be extracted directly from observations, rather than 
from simulations. 
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Figure 10: A graphical depiction of the regions in the {k,z) plane corresponding to the validity of one- 
and two-loop EFT, the paths along which the integral for the lensing potential power spectrum Cf is 
evaluated, and the scales and redshifts that contribute 95% of the value of Cf for selected £ values. We 
find that one-loop EFT can be used to calculate Cf with less than 5% theoretical error for £ < 600, while 
the corresponding limit for two-loop EFT is £ < 1000. See the main text for further discussion. 


over Pff) times a geometric factor: 

= . (32) 

where x is fh® comoving distance with respect to a chosen observer, running from zero to x*) 
the distance to the emitting source (for CMB lensing, the source is located at the last scattering 
surface). The derivation of this formula, as well as a thorough discussion of gravitational lensing, 
can be found in [24]. The EFT prediction for Cf can be calculated by using Peft-i- 1 oop or Peft- 2 - 1 oop 
for Ps in Eq. (30), and then inserting (30) into (32). 

To investigate what ranges of £ can be accurately described using the EFT at one or two loops, 
we plot several sets of curves in Fig. 10 for lensing of the CMB. The solid lines trace the paths in 
the {k, z) plane along which the integral (32) is evaluated for various values of i, while the black 
short-dashed line separates each path into the region contributing 95% (above the black line) or 
5% (below the black line) of the value of the integraT®. This shows, for example, that Cf^^QQ is 
dominated by the region of the integral corresponding to /c < 0.25 h Mpc“^ and 2 ; > 0.55, while, as 
expected, higher-fc modes contribute more strongly at higher £. The grey long-dashed lines show 
where the one- and two-loop EFT predictions for the matter power spectrum fail, as found in 
Sec. 3.4 with WMAP7 cosmological parameters. For each £ curve, if the 95% region overlaps with 

^®We use Halofit [23] for this calculation, since it provides a rough indication of the size of nonlinearities that 
will affect the boundary between 95% and 5% contributions to the integral. 


24 







Figure 11: As Fig. 10, but for the tensing potential power spectrum corresponding to sources at z = 2 
(top panel) and z = 1 (bottom panel). The two-loop EFT prediction for Cf has less than 5% theoretical 
error for i < 600 for sources at z = 2 and I < 350 for z = 1. 


the region where the one- or two-loop EFT fails, then the theoretical error in the EFT calculation 
of cf will surpass 5%, which we take as our threshold for a useful calculation^®. This occurs at 
£ ~ 600 for one loop and £ ~ 1000 for two loops. 

In Fig. 11 we repeat Fig. 10 but for sources at redshift z = 2 and z = 1 corresponding roughly 
to the reach of ongoing lensing surveys. For photons originating at these redshifts (as opposed 
to the last scattering surface), the effects of late-time nonlinearities are more important. We see 
that the EFTofLSS allows us to predict Cf to 5% accuracy up to £ ~ 600 for sources ai z = 2 
and £ ~ 350 for z = 1. This is a large improvement with respect to former analytical techniques. 

In Fig. 12, we plot the ratio of the two-loop EFT and SPT predictions for Cf over the pre¬ 
dictions from linear theory along with grey bands indicating the theoretical error on the EFT 

Observational uncertainties on future measurements of Cf will precisely determine the required theoretical 
precision. For now, we choose 5% as an optimistic forecast of these future requirements. As presented in recent 
publications, current measurements of Cf bandpowers have uncertainties that are well above this level [25, 26, 27, 
28] , although these uncertainties could be reduced by changing how the measurements of individual multipoles are 
binned. 

^°In computing the two-loop EFT prediction for Cf , it is important to recall that Feft-2-1oop is written as an 
expansion in powers of k/kjsih times Pn, and at high k these terms will cause the prediction to deviate much more 
strongly from Pnl than Pn does. Such a deviation is guaranteed regardless of the loop order, since as k ^ /cnl the 
entire perturbative expansion breaks down completely. To avoid using the EFT prediction in a regime where it does 
not apply, for the EFT curves in Fig. 12, we use Peft-2-1oop for Ps when k < kfjfz), and Pn when k > k^Jfz). 
This introduces a small error in the calculation of Cf , which should be no more than 5%, since we have already 
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Figure 12: The ratio of the two-loop EFT and SPT predictions for Cf to the prediction from linear 
theory. The grey bands show the theoretical error on the EFT prediction, based on the contribution to the 
cf integral from scales beyond the reach of the EFT. It is likely that future data will become sensitive to 
the nonlinearities in the lensing potential that are encoded in the differences between these curves (which 
may also be accessible to current measurements of CMB lensing, after an appropriate choice of binning 
scheme). We furthermore see that, for I > 200, SPT at two loops is as bad an approximation to the 
EFTofLSS (which we showed is the correct result) as linear theory. 


prediction, based on the contribution to the lensing integral from scales beyond the EFT’s reach. 
As presented in recent publications, current data are unable to distinguish between these curves, 
since, with the chosen binning schemes, the errorbars on Cf are 0(10%) [25, 26, 27] or slightly 
smaller [28]. For example, the 2015 Planck data [27] measure the overall amplitude of the lensing- 
potential power spectrum in the range 40 < £ < 400, to a precision of about 2.5%. Fig. 12 
shows that this might allow us at most a very marginal measurement of gravitational nonlineari¬ 
ties Future data are expected to be precise enough to detect these effects with high signihcance, 
making the EFTofLSS an important tool for using lensing to tighten constraints on cosmological 
parameters and search for new physics. 

Thus far, we have presented results for dark matter only, so an important question is whether 
the effects of baryons could become important at the scales we are concerned with. At lowest 
order, these effects will manifest in the power spectrum as a shift of ~15% in the value of 

in the tree-level counterterm (see [14] for a detailed discussion). To quantify the impact on 
the lensing potential power spectrum, we calculate the change in Cf resulting from a 15% shift 
in c^(i)(^) at each z, and compare to the full two-loop Cf calculation. Fig. 13 shows that, in the 
range of multipoles where the EFT prediction has better than 5% theoretical error, and for the 
sources we consider, baryonic physics never alters the lensing potential power spectrum by more 

established that for I < 1000 the region k > contributes no more than 5% of the value of the integral 

corresponding to lensing of the CMB. Using Peft instead of Pu in this regime would probably make the error 
induced by the contribution from high wavenumbers more important. 

Notice that one could design a procedure that would allow us to measure an integral over redshifts of the speed 
of sound, which would allow us to extract this information directly from data, without relying on simulations. 
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Figure 13: The effect of baryonic physics on the two-loop EFT prediction for Cf, as approximately 
captured by a 15% shift in in the tree-level counterterm Pfrel [H]- The red arrows indicate the 

ranges in which each EFT prediction has less than 5% theoretical error, as estimated from Figs. 10 and 11. 
We find that baryonic effects on Cf never exceed r^l% in these ranges and for these sources. 


than ~1%. 

In summary, all of this is a promising indication that the EFT prediction could be applied to 
upcoming observational data, an exciting prospect that we plan to explore in future work. 

5 Discussion 

In this paper we have studied the prediction of the EFTofLSS as applied to dark matter at 
different redshifts. This has allowed us to explore to which high wavenumber the EFTofLSS can 
push itself as we increase redshift. We have found that the reach of the EFT in the UV grows 
with redshift, largely outperforming SPT and all former analytic techniques, all of which, when 
correctly implemented, have the same UV reach as SPT. At two loops with IR-resummation, the 
EFTofLSS agrees within approximately 1% with matter power spectra measured from simulations 
up to fc ~ 0.6hMpc“^ at z = 0. At redshifts z > 1, the agreement is within approximately 
2% up to at least k ~ IhMpc”^, and also within approximately 5% up to A: ~ 1.2hMpc“^ at 
z = 1 and k ~ 2.3hMpc“^ at z = 3. Uncertainties in the htted value of and the range of 
reasonable choices for fcren can change where the predictions fail; these uncertainties are largest at 
low redshift, where the calculation is particularly sensitive to accidental cancellations that may 
occur between different terms that enter at two loops, and become milder with increasing redshift. 

As we go to higher redshift, the slope of the power spectrum near the nonlinear scale becomes 
closer to —3. This makes the effect of loops more important relative to the effect of the countert¬ 
erms, so that the gain of the EFT with respect to SPT decreases at higher redshifts. In fact, in 
the EFTofLSS it is possible to estimate the size of the higher-order contributions that have not 
been included, so that one can forecast the scale at which a calculation at a given order fails. 
After pointing out that approximating the universe as scaling is only marginally correct, we have 
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developed a more accurate way of estimating the size of higher order terms, with the results that 
the UV reach of our two-loop calculation (along with the associated uncertainty in that reach) 
agrees with these estimates. 

We have studied how the parameters of the EFT depend on time, Ending that, in the local 
in time approximation, the term can be well described by a function close to a power-law 
that slightly runs with z, and whose generic features can be anticipated by using an approximate 
scaling symmetry of the universe, and that ultimately requires a single parameter to be fit. In 
summary, only two coupling constants—at redshift z = 0, and the slope of the power law, 
controlled by a parameter f3 —are sufficient to extend the predictions of the power spectrum in 
the EFTofLSS to all redshifts far beyond the reach of SPT. 

These power spectrum results allow us to perform analytical calculations related to gravi¬ 
tational lensing, such as the lensing potential power spectrum presented in this work. We have 
found that the EFTofLSS at two-loop order is able to predict Cf corresponding to lensing of CMB 
photons with better than 5% accuracy for i < 1000. For photons emitted by galaxies at z = 2 or 
z = 1, similar accuracy is achieved by the EFT prediction for £ < 600 and 350, respectively. This 
opens up an incredibly interesting window to compare the EFTofLSS to currently available data 
as well as to data that will be available in the very near future. We plan to perform a detailed 
study in an upcoming paper. 

We stress that the presence of coupling constants in the EFTofLSS is not a catastrophe for 
the theory. These parameters can be measured directly in observations (without resorting to 
simulations), or can be measured in small (and therefore fast) numerical simulations It would 
be interesting to explore how much cosmological information is lost if these coupling constants 
were to be measured directly in observations. We guess not much, given that in the EFTofLSS the 
same coupling constant predicts several observables: fixing the speed of sound parameter using the 
matter power spectrum allows us to predict the matter bispectrum [10] and the momentum power 
spectrum [9] at one loop up to about k ~ 0.3hMpc“^ without any remaining free parameter. We 
leave a more detailed study of these aspects to future work. As we keep exploring the predictions 
of the EFTofLSS for dark matter, biased tracers, lensing, and so on, the repeated successes of this 
novel analytic approach seem to indicate that we might be able to extract a much greater amount 
of cosmological information from large scale structures than previously believed. 
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Appendix 

A Derivation of Tree-Level Counterterms 


We will now present the derivation of the tree-level counterterms involved in one- and two-loop 
EFT predictions, since one of these terms was erroneously omitted from prior work [3, 9]. For 
convenience, we will repeat some of the equations from Sec. 2. 

We begin with the Eulerian-space EFTofLSS equations of motion, with the stress tensor 
given by Eqs. (4) and (7): 


oH daS{k, a) + 6{k, a) = — a{q, k — q)S{k — q, a)6{q, a) , 


(33) 


oH daO{k, a) + l-i6{k^ a) -|- -'H^(a)f2m(a)(5(fc, a) = — /3{q, k — q)6{k — q, a)6{q, a) 

^ iqr 

da' 




(34) 


where e is a parameter inserted to organize the power of K{a,a') appearing in the perturbative 
solutions we will obtain. To solve these equations, we make use of two assumptions: 

• ffm(a) ~ /(a)^ and 

• Cn{a) = where 


Cn{a) 


da' 

a'H{a') 


K{a, a') 


DijaT 

Di{a)^ 


(35) 


and ^ is a constant that we are free to set to a convenient value, since it simply rescales 
each c„. 


Solutions to (33) and (34) can also be found without making use of either of these assumptions 
(e.g. [16]), but we use them here because they greatly simplify the algebra involved, and we 
have checked that the exact implementation of the time-dependence does not make a relevant 
difference [16]. Under these assumptions, we can use the following solutions: 


S(a, k) = X^[B,(a)]"+"<i'“>(fc) , 

n=l n=l n=l 

f OO CO oo 

«(o, k) = -H{a)!(a) \ j;iD,(o)|"e<”l(fc) + e5;ifl.(a)|”+<«'”)(fe) + ^|B,(o)]”+^«e<”>(A;) 


n=l 


n=l 


n=l 


(36) 


Since in this Appendix we are only interested in tree-level counterterms, we will only need the 
n = 1 terms in each sum above. Plugging these back into (33) and (34) and collecting terms of 


29 





order e^, we find: 


+ = 0 , 


- (^ + c) e">(fc) + p"’(fc) = 4"c,5i'‘>(fe). 


( 37 ) 


Solving for S^^\k) and conveniently setting ^ = ('((' + 5/2) to cancel the resnlting denominator, 
we find 

^(i)(fc) = -ciPsW(k) . (38) 

In [3], Cl was identified with (27r)c^^j^^//cNL5 since the calcnlation was only performed at a single 
redshift. In this paper, it will be convenient to absorb a factor of [Di(z)]‘’ into a time-dependent 
fnnction c^(i)(^), coinciding with the convention from [3] at z = 0. Therefore, the connterterm 

corresponding to 2 (which enters the one-loop EFT prediction) is written as 

z){k, z) = -2{2Tr)ci^^{z)[Di{z)]‘^^Pn{k) , (39) 

%L 

while the connterterm corresponding to takes the form 

2 ) = (.2wn4m(z)riD,(z)f-Ppn(k) . ( 40 ) 


The second connterterm that scales like {k/k^i^)^Pii arises from correlating with To 
End the expression for we follow the same steps that lead to Eq. (37), but collect terms of 
order instead of e^; this gives 

(1 + 2C)S^^\k) - 0«(fc) = 0, - Q + 2C) 0W(fc) + ^<5W(fc) = k%+^^S^^\k) . (41) 

Solving for m^)(fc) and using our previous choice for we hnd 

Sf'\k) = (k)^ ( 42 ) 


At this point, we will assume locality in time, for which = ci [3]. Then, taking 2(5*^^^^^^) 
gives us our second connterterm: 

PP(k,z) = {2n)^:^^y.^,-,(z)nDdz)?Ppu{k) . ( 43 ) 

+ 4 ) %L 

The (/-dependent factor in Ptree’^^ is a slowly-varying function of ( for 1 < (/ < 5, so deviations 
from our initial choice of (/ = 3 will have only a weak effect on this term. 
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Figure 14: Exact and approximate versions of /Pn at z = 0, where P^^^op^{k,0) = P 2 -ioop{k) — 

2(27r)Cg^2)(0)(^/^NL)^Fii(A:). The estimate is plotted using Eq. (22) with c^(;^)(0) = 0, a ~ 3/10 and 
ft = 1. The two curves match to within a factor of ^2 up to k ^ 0.7 /iMpc“^; an even better match can 
he achieved by letting a or fd float freely instead of fixing their values to those from Sec. 3.2. 

B Notes on Estimates 

B.l Estimating Finite Parts of SPT Loop Corrections 

In this short section, we simply highlight the fact that the method of estimating loop corrections 
from Sec. 3.2 also applies to the hnite parts of SPT terms alone. At two loops, this corresponds 
to the combination 


dilpA. 0) = J’2-io.p(fc) - 2(2^)c2p|(0)(4/4NL)"ni(i) , (44) 

where the c^^ 2 ) term removes the UV contribution to P 2 -I 00 P (see Sec. 3.3). Fig. 14 compares 

the exact calculation of with the approximation from Eq. (22), with c^(i)(0) = 0 and a 

and ft hxed to the values used in Fig. 4 (3/10 and 1, respectively). The exact and approximate 
calculations match to within a factor of ~2 up to /c ~ 0.7 hMpc“^, with the discrepancy due 
entirely to the normalization of the estimate, set by the combination /? x a in this case. Letting 
ft or a float away from the values used in the plot enables an even better agreement. Since the 
SPT contribution is a part of the EFT one, it is important to be able to estimate that as well. 

B.2 Comparisons of Estimates and Exact Calculation of i^ 2 q*oop^ 

For the interested reader, in this section we provide comparisons between the estimates for T’ 2 noo^p^ 
from Sec. 3.2 and the exact calculations from 2 ; = 1 to 4. As with Fig. 4, in these plots we have 
hxed to its best-ht value at each redshift, and have hxed a to 3/10 plus any small shifts that 
are required to remove zero crossings from the estimated curves. One can see from Fig. 15 that 
the estimates provide a good match to the exact calculations at all redshifts, though one should 
keep in mind the presence of the same uncertainties that we describe in Sec. 3.2. 
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Figure 15: Comparisons of estimates and exact calculations of P 2 - 1 I 0 P for z = 1 to A. The estimates 
are computed as described in Sec. 3.2. 


C Additional Discussion of 


As mentioned in Sec. 3.3, it is possible to calculate analytically in a scaling universe by taking 
the limit of P 2 -I 00 P in which the external momentum k is much smaller than the internal momenta 
q and p. Explicitly, if Pii{k) oc fc”, then this limit of the integral will be dominated by the term 
containing the kernel, which scales like in the k q limit. Therefore, the UV limit of 

F 2 -ioop(fc) looks like 



-r^Pnik) = 


k^ 


Piiik) 



(45) 
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Figure 16: We plot the full P 2 -I 00 P calculation in ACDM (solid), a c^( 2 )(^/^nl)^Fii term with c^^ 2 ) M 
at kren = 0.05/iMpc“^ (dotted), and the sum of two terms (dashed): a Pii term with c^^ 2 ) 

fit at kren = 0.3/iMpc“^, and a log{k/kren){k/kt^i^yPii term. The dashed curve ean be made to agree 
with P 2 -I 00 P at low k if the coefficient of the log term is set to C\og —0.55. This implies that we cannot 
exclude the possibility that P 2 -I 00 P contains a contribution that depends logarithmically on k, and therefore 
c^( 2 ) may depend on kren- 


where C* is a constant that can be calculated exactly by evaluating the q and p integrals, and 
e is an order one number. (The 0{k/k-i^if) factor comes from the lower limit of integration, but 
evaluating (45) with k /cnl ensures that this factor will be highly suppressed, allowing for an ac¬ 
curate determination of C.) Setting c^|. 2 ) oc —(after restoring the proper conventions) will 
precisely cancel this contribution at all values of k. Alternatively, one could £x c^^ 2 ) by matching 
to the full calculation of P 2 -ioop{k) at some small but hnite value of k, where the (fc/fcNL)^A’ii(fc) 
term in P 2 -ioop{k) is dominant. Note that this entire discussion is taking place at 2 ; = 0, but the 
time-dependence of cl^ 2 ) ''^^1 be hxed to the time-dependence of P 2 -I 00 P, so there is no need to 
distinguish between different times. 

The situation changes slightly at certain values of n, however. A case of interest is n = —2, 
which is quite close to the slope of the power spectrum on quasi-linear scales in our universe 
for which depends logarithmically on A: 


p(UV,n=-2) 

■^2-loop 


(fc) 


C'fin T l^log lo§ 




Pii{k) . 


(46) 


In order for c^^ 2 ) cancel this A-dependence, it will gain a dependence on a renormalization 
scale kren- c^( 2 )(^ren) [—Cfin + Clog log(A//cren)]• The Same phenomenon occurs when loops are 
renormalized in quantum held theory, where the renormalization scale is often denoted by p and 
Clog is the analogue of a /^-function. If we were to attempt to determine c^( 2 ) by matching to P 2 -I 00 P 
at low k, the value we would hnd would depend on the scale fcren at which we perform the matching, 
with the difference between matching at two scales fcren.i and fcren ,2 given by Cog log(fcren,l/fcren, 2 )• 

^^For this discussion, the relevant scales are the ones close to where we choose the renormalization scale. 


33 















z 


Figure 17: As Fig. 7, but for WMAP3 instead of WMAP7 cosmological parameters. The red curve 
shows the fitting function (25) with /3 = —1.20, which exhibits excellent agreement with the values of 
c^(i)( 2 ;) obtained separately at each redshift. 

In ACDM, the effective slope rieff (fc) of the linear power spectrum varies strongly with k, so that 
the UV part of cannot be determined as cleanly as in the scaling case. Since Uefr(fc) passes 

through —2 around k ~ 0.25hMpc“^ (see Fig. 1), we cannot exclude the possibility that 
contains some kind of logarithmic dependence on k around this scale, implying that 0 ^( 2 ) would 
depend on the scale at which it is fixed. To show this concretely, Fig. 16 displays the full P 2 -I 00 P 
curve, along with two other calculations. The dotted curve corresponds to c^(' 2 )(^/^nl)^-Pii where 
c^(^ 2 ) has been matched to P 2 -I 00 P at ^ren,! = 0.05hMpc“^. The dashed curve shows the sum of 
'^s( 2 )(^/^nl)^-Pii with c^^ 2 ) fit using the procedure from Sec. 3.3 (with an effective renormalization 
scale of /cren ,2 = 0.3hMpc“^), and a log term corresponding to renormalization at /cren, 2 - The 
fact that both curves match equally well with P 2 -I 00 P at low k shows that we cannot distinguish 
whether or not there is a log term in P 2 -I 00 P; and therefore we cannot exclude the possibility that 
c^( 2 ) depends on fcren- 

Therefore, to be consistent, we must use (roughly) the same scale to determine both and 
c^( 2 ). We cannot use data at very low k to determine due to the fact that at low k the 
contribution of the term becomes smaller than the numerical, and cosmic variance, error in 
the simulations. This implies that we should not match to P 2 -I 00 P at low k, but should instead 
use information at higher scales to determine c^( 2 ) Tfi® precise procedure we use, based on 
matching at a certain z-dependent fcren, is described in Sec. 3.3. The fact that fcren depends on 
z also explains why the time-dependence of c^^ 2 ) i® aot simply given by the time-dependence of 

P 2-loop • 
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Figure 18: As Fig. 8, but for WMAP3 instead of WMAP7 cosmological parameters. The reach of the 
EFT predictions is comparable what we find in Sec. 3.4- 

D Consistency Check using Different Cosmology 

In this appendix, we apply the procedures in Secs. 3.3 and 3.4 to the nonlinear power spectrum 
generated by Coyote using WMAP3 cosmological parameters: h = 0.73, = 0.238, fib = 0.0418, 

fiA = 1 —fim, = 0.951, and ag = 0.74. Fig. 17 shows the values of determined separately 

at each redshift, along with the fitting formula (25) with the best-fit value f3 = —1.20, while Fig. 18 
compares the EFT predictions with the nonlinear spectra at different redshifts. Overall, we find 
results that are comparable to what we found for a WMAP7 cosmology, both for the reach of 
the EFT prediction at various redshifts and for the validity of Eq. (25) for describing the time- 
dependence of 

^^Alternatively, one could measure Cg^ 2 ) s-t low fc’s and then compute the running by numerically evaluating the 
two-loop integral, an approach that seems numerically quite challenging. 
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E Effect of Non-Locality in Time 

We have performed our main study using the assumption that the short modes that have been 
“integrated out” of the EFTofLSS affect the long modes in a way that is purely local in time, an 
assumption that may not be strictly correct in the real universe [3, 7, 12], This assumption affects 
the two-loop matter power spectrum through the terms and ■ The first of these can 

be written in terms of three separate functions Pn{k) and three coefficients c„ (see [3] for more 
details): 

= Ci(z)Pi(k) + c,(z)P,(k) + cz(z)p,(k) . (47) 

while the second can be written as 

Piref\k,z) = ^^^^Ci{z)ci+^{z)[Di{z)]‘^k^Pn{k) . (48) 

The hrst coefficient, Ci, is dehned by Ci{z) = (27r)c^(-^^(2:)//cNL) while c„ = if Eq. (7) is 

used for the kernel K{a, a') appearing in the Euler equation. The parameter p is a simple way to 
parametrize the severity of the non-locality in time: a lower value for p corresponds to a broader 
kernel, while p —)■ oo corresponds to the local limit, effectively turning the kernel into a delta 
function at a' = a. 

In Fig. 19, we repeat the procedures in Secs. 3.3 and 3.4 using different values of p in the 
two-loop EFT prediction, re-fitting in each case. The effect of changing p is most prominent 
a.t z = 0, when the greatest amount of time evolution has taken place. For the most nonlocal 
case we consider, p = 1, there is a ~2% shift in the power spectrum for k > 0.2hMpc“^, with 
a corresponding 6% shift (upward) in the value of c^(i)(0). For higher, and maybe more reason¬ 
able [12], values of p, the effect on the power spectrum is much less, especially when compared 
to the error induced by the uncertainty in (see Fig. 8). The non-locality has even less of an 
effect at higher redshift, and, importantly, cannot improve the reach of the prediction into the UV 
at any redshift. 

Instead of using the relationship between the c„ coefficients implied by Eq. (7), we can also 
allow for a completely general type of non-locality by treating each c„ as independent. The 
two-loop EFT prediction then has four free parameters at each redshift, which we can write as 
c^(i)(^) and rn{z) = Cn{z)/ci{z) for n = 2, 3,1 -|- C- We determine the four parameters at each 
redshift through a least-y^ £t to the nonlinear power spectrum, subject to the constraints that 
ki+cl — I^31 — I^ 2 1 < 1 and cl^i){z) is within 30% of the value obtained from the time-local fit. 
The first constraint comes from the fact c„ is equal to an integral over the power of the growth 
factor Di{z), and 0 < Zli(z) < 1, so higher powers of Di{z) in the integrand will act to suppress 
the absolute value of the integral. The second constraint comes from the fact that can be 
determined from the one-loop calculation, and the value for when adding time non-locality 
at two loops, should not be drastically different from the one-loop value. 

At each redshift, we perform this fit between fc m^n = 0.15 hMpc“^ and a ffoating fcmax, taking 
the overall best fit to be one that gives a fiat region in PEFT- 2 -ioop/-Pcoyote that reaches farthest 
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Figure 19: Two-loop EFT predictions for the WMAP7 power spectrum, parametrizing the severity of the 
non-locality in time of the counterterm with a single parameter p. The speed of sound parameter 

c^(i) is fit separately at each redshift and p value using the procedure described in Sec. 3.f. Taking the 
uncertainties on the nonlinear data into account, we find that introducing non-locality in time does not 
improve the agreement of the prediction with the data at any redshift. 


into the UV. (For the hts at z = {0,1,2}, the fcmax values corresponding to these best hts are 
(0.6,1.4,1.9} hMpc”^.) The results are shown in Fig. 20. We hnd only minor differences between 
these general predictions and those of the time-local approximation. In summary. Figs. 19 and 20 
demonstrate that the matter power spectrum at z < 2 (and also at 2 ; < 4, although we do not 
plot redshifts greater than 2) does not provide any compelling evidence for non-locality in time, 
although such evidence may emerge from consideration of other observables. This is an interesting 
topic for future studies. 

We finally point out the following. Since, from comparison to data, we find no need to abandon 
the time-local approximation, the EFTofLSS at two-loop order has only two free parameters for all 
redshifts: at 2 ; = 0 and its ‘ 2 :-slope’ as described. In the non-local in time case, additional 

parameters are present, but it is incorrect to consider these additional parameters as having been 
fitted to obtain the results presented in the main text of this paper. This is because the local-in¬ 
time case constitutes a special, highly symmetric, point in the theory space of the EFTofLSS. As 
is standard, when theories are enhanced by additional symmetries or properties, the number of 
free parameters is reduced. 
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Figure 20: Two-loop EFT predictions for the WMAP7 power spectrum, allowing for a generalized form 
of non-locality in time in which the four Cn coefficients appearing in and Plref'^ (VO 

and (48)) are treated as independent, and are fit as described in the main text. This results in predictions 
that are not significantly different from the time-local case, motivating the adoption of the time-local 
approximation in practical computations. 
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